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Properties of spin-1 Bose gases with ferromagnetic interaction in the presence of a nonzero mag- 
netic field are studied. The equation of state and thermodynamic quantities are worked out with 
the help of a mean-field approximation. The phase diagram besides Bose-Einstein condensation 
contains a first order transition where two values of the magnetization coexist. The dynamics is 
investigated with the help of the Random Phase Approximation. The soft mode corresponding to 
the critical point of the magnetic phase transition is found to behave like in conventional theory. 
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I. INTRODUCTION 



Atomic Bose gases consisting of atoms with spin 1 can 
exhibit a spinor condensate, which has attracted consid- 
erable attention in recent years 00ESHIE0,H&[I3, 
0,0. Such gases have interesting properties also before 
the Bose-Einstein condensation (BEC) sets in. This is 
particularly true when the coupling between the spin de- 
grees of freedom prefers ferromagnetic ordering. Such a 
system can be realized by the gas of 87 Rb atoms. A 87 Rb 
atom has a nuclear spin of j = 3/2 and an electron spin 
of s = 1/2, therefore its net spin can be F = 1 or F = 2. 
Since the energy of the F = 1 multiplet is smaller than 
that of the F = 2, the population in the F = 2 multiplet 
dwindles at sufficiently low temperatures resulting in a 
gas of purely spin-1 (3-component) bosons 0,0- 

Az zero external magnetic field it has been found that 
when lowering the temperature a transition to ferromag- 
netic state can occur first which can be first or second 
order as well depending on the strength of the spin-spin 
interaction The transition to the Bose condensed 

phase can also be first or second order according to mean- 
field theory results. Note that the tendency towards fer- 
romagnetic ordering is present in the Bose gas formed 
by atoms with nonzero spin even in the case when only 
a spin independent interaction acts between the atoms 
(or the gas is ideal) |14L Il5l 1 161 Il7| . Effects of an ex- 
ternal magnetic field has been studied for the free Bose 
gas in Ref. 0. In the present paper a gas of spin-1 
bosons is investigated in the presence of an interaction of 
ferromagnetic type and a nonzero magnetic field. 

The system is supposed to be translationally invariant 
with a homogeneous magnetic field pointing to the z- 



directionand. The Hamiltonian takes the following form: 
"H = ^ [( ek - ^ Srs ~ (F z ) rs <4(k)a s (k) 

k 

a^(ki)4(k 2 )K;' s ' as (k3) as ,(k 4 ), (1) 



k l+ k 2= k 3+ k 4 

where a£(k) and a r (k) create and destroy one-particle 
plane wave states with momentum k and spin projection 
r. The spin index r refers to the eigenvalue of the z- 
component of the spin operator and can take values from 
+, 0, — . In this basis the spin operators are given by: 
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In Eq. (JTJ ek = H 2 k 2 /2M refers to the kinetic energy of 
an atom (M is the mass of an atom), /1 to the chemical 
potential, g to the gyromagnetic ratio, /ib to the Bohr 
magneton, B to the modulus of the homogeneous mag- 
netic field, V is the volume of the system and V^ s s the 
Fourier transform of the two particle interaction poten- 
tial, which for the low temperature, dilute gas can be 
modeled by the momentum independent s-wave scatter- 
ing amplitude given for spin-1 bosons by [Tsl Il9(: 



V, 



iS rs S r ' s f -f~ c s (F) rs (F) r '_ 



(3) 
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with parameters: 



Airh 2 a + 2a 2 
~M 3 : 
AttH 2 a 2 — a 
~M 3 ' 



(4a) 
(4b) 



The parameters ao > and 0,2 > are the scattering 
length in the total hyperfine spin channel zero and two, 
respectively. Note, that c s < for the gas of 87 Rb atoms 
|20|. For such a system it is energetically favorable to 
align the spins along one direction, i.e. the system has 
a ferromagnetic coupling 0, Q • In this paper we assume 
such a system. 

The outline of the paper is as follows. In section 
ITT1 the Hartree equation of state of the dilute and low- 
temperature, interacting spin-1 Bose gas is given in the 
presence of a nonzero magnetic field. In section ITTT1 the 
magnetic susceptibility of the system is investigated in 
the uncondensed phase. Section Hvl is devoted to the de- 
termination of correlation functions of generalized den- 
sity operators, which are used in section[V]to express the 
linear response functions of the system. With the help of 
the response functions further static properties and the 
spin dynamics of the system are also investigated. In 
Section IVII the results are summarized. 

II. EQUATION OF STATE 

For the spin-1 Bose gas the Hartree approximation [2l| 
yields a plausible set of equations, which form the equa- 
tion of state of the system. The first equation in the set 
expresses the total density as the sum of the density of 
the condensate and that of the different spin projections 
of the non-condensate and reads as: 
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n 



(5) 



where n = N/V the total density of atoms in the gas, 
?io = Nq/V the density of the Bose- Einstein condensed 
atoms, and 



k 



l k,r' 



(6a) 



is the density of the noncondensed atoms with spin pro- 
jection r, where 



n k,r — 



(6b) 



£k,r = e k 



/i + c n n + r(mc s — g^B). (6c) 



In the last equation m is the magnetization density of 
the system [see below] and (3 — l/k^T is the inverse 
temperature. The term rmc s can be interpreted as the 
energy shift of the internal state r arising from a "molec- 
ular field" . The formula is quite plausible and can also 
be derived from a self-consistent Hartree approximation 



in the Green's function technique. The next equation of 
the set is the expression of the magnetization density m, 
and it takes the following form: 



m = no + ti 1 — n_ 



(7) 



since the condensate spinor points to the z-direction. In 
the case, when uq = 0, i.e. the system is not Bose- 
Einstein condensed Eqs. © and (JJJ) together with the 
expressions (jHJl form a closed set of equations, and in 
the knowledge of the temperature, magnetic field and 
particle density can be solved to the chemical poten- 
tial and the magnetization density. However, with low- 
ering the temperature, one arrives at a point, when 
the lowest bound of the energy expression will be zero: 
e i + = 0, i.e. the chemical potential will reach the 
value of n = c n n + c s m — g^B, which means that 
n' v , = [exp(/3ek) — 1] — x , so the system undergoes Bose- 
Einstein condensation. This remains true at lower tem- 
peratures, i.e. 



= ?io [ — /i + c n n + c s m — g/isB] 



(8) 



holds. The multiplicative factor no is used to make the 
equation valid for the high temperature phases as well, 
where no = 0. In conclusion Eqs. JSJ), and (JSJ to- 
gether with the expressions © form a closed set of equa- 
tions for all possible temperature values and is considered 
as the equation of state of the spin-1 Bose gas with ferro- 
magnetic interaction in the presence of a magnetic field. 
Our equation of state compares for zero magnetic field 
with that obtained by Gu and Klemm |13| in a com- 
pletely different approach, while for c s = it becomes 
similar to that of the free Bose gas treated in a nonzero 
magnetic field by Simkin and Cohen |l7j . 

For the solution of the equation of state let us restrict 
ourselves for the fixed particle density, magnetic field and 
temperature case, i.e. Eqs. ©, Q and © are solved 
for (ji, m, no) with (n, B, T) fixed. It is important to note 
that with introducing // = /i — c„n the effect of c n n can 
be incorporated to the chemical potential for fixed par- 
ticle density. Detailed investigations show that the mag- 
netic separation found is accompanied by the separation 
of density. The density difference, however between the 
phases is smaller by a factor of order |c s |/c n than the 
density itself. To concentrate to the magnetic properties 
alone we disregard the phase separation in the density. 
The phase diagram including phase separation will be 
published elsewhere [22|. 

The BEC transition temperature of the homogeneous 
spin-1 Bose gas is given by 



To 



2nh 2 



3C( 



(9) 



and by taki ng Tp t o be, e.g. 200 nK (as is typical in 
experiments |2J,|2J|) the particle density is given by fun- 
damental constants (apart from the mass of the atoms in 
the gas). For the case of 87 Rb atoms, with T = 200 nK 
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FIG. 1: Possible phase diagrams of the system: a) for e s > e» , b) for ei 2 ^ < e s < el , c) for ei 3 ' < e 3 < ei 2 \ d) for e s < e. 
The dashed line symbolizes a first-order phase transition, while the continuous line refers to a continuous phase transition. 



(3) 



the particle density by Eq. J^Jl is n — 10 20 m~ 3 . It is 
convenient to introduce the following dimensionless pa- 
rameter: 



n\c s \ 



(10) 



which can be interpreted as a sort of mean-field energy 
in units of fee To- 

By solving the equation of state the domain of e s can 
be divided into four parts according to the character 
of the occurring magnetic and BEC phase transitions. 
In the first region, when e s is greater than a certain 



value {e s > e 



(i) 



1.22) there exist a paramagnetic- 



ferromagnetic transition for B — at T — Tm above the 



BEC transition Xbec • Such a phase diagram is depicted 
in Fig. da) f° r £ s = 1-4. Both transitions are continu- 
ous for these values of e s ; the magnetic transition exist 
only in the absence of magnetic field and its critical tem- 
perature is higher than the BEC transition temperature 

If the value of e s is smaller than ei 1 ^ the magnetic tran- 
sition at B = becomes a first order one. Therefore the 
point (e s = ei 1) ,S = 0) can be understood as a tricritical 
point for the magnetic transition, since it separates the 
region of first-order and continuous transitions. The first 
order transition survives even in the presence of a small 
magnetic field leading to the coexistence of two phases 
with different magnetizations. The transition, called the 



4 




is shown in Fig. \I\d) for e s = 0.6. 

After discussing the above four possible situations (ac- 
cording to the value of e s ), it is worthwhile to plot the 
full phase diagram of the system, where e s is also con- 
sidered as a variable besides T and B. It can be seen in 
Fig. [3 The above four figures can be obtained as sec- 
tions of Fig. [21 made with a fixed e s . It is important to 
note that inside the simply- hatched region on the (e s , B) 
plane of Fig. [3] the magnetic transition is of first order 
and its transition temperature is higher than that of the 
continuous BEC (of the same e s and B value) . 



III. MAGNETIC SUSCEPTIBILITY IN THE 



UNCONDENSED PHASE 



FIG. 2: The phase diagram of the spin-1 Bose gas in the 
space of (e s ,B,T). The temperature is made dimensionless 
with the help of To, while the magnetic field is given in units 
of /cbTo/<?Mb. The solid surface indicates a continuous BEC 
transition, while in the cross-hatched region BEC is of first 
order. The simply hatched region, outside the region of first 
order BEC, refers to the first order magnetic transition and 
is lifted slightly from the surface of continuous BEC. At the 
heavy curve the magnetic transition is continuous, and its 
transition temperature is above that of the BEC. The sur- 
faces of first order transitions and the curve of the continuous 
magnetic transition is projected to the (e s , B) plane for en- 
lightenment. The relevant e s values are also indicated [com- 
pare with Figs. Qa)-d)]. 



magnetic transition (MT) in the following, ends in a mag- 
netic critical point (MCP), specified by the critical value 
of the temperature T c and that of the magnetic field B c 
for fixed e s . Note that a first order phase transition 
due to an external field coupled to the order parame- 
ter was observed in the ferroelectric Ba Ti O3 already in 
1953[2! H3- If at the same time e s > ei 2) w 0.78, BEC 
remains a continuous phase transition and the phase di- 
agram of this region looks like plotted in Fig. Qb) for 
e s = 0.9. 

If e s becomes smaller than , BEC becomes also of 
first order at B — 0. However for g /ie-B / 'fceTb ^S> 1 the 
system behaves like the gas of noninteracting scalar par- 
ticles with a continuous BEC phase transition at 3 2 / 3 Tb, 
therefore there must be a critical magnetic field value, 
B^\ which is also a tricritical point, but in the sense 
of the BEC transition, which separates the first order 
BEC in small magnetic field from the continuous BEC in 
strong magnetic field. Such a phase diagram is plotted in 

Fig. |l]c) for e s = 0.72. As long as e s > e { s 3) « 0.62 there 

(2) 

exist another critical magnetic field value (Be ) above 
which there is no magnetic transition (magnetization be- 
haves analytically). 

(3) 

For e s < e s the two critical magnetic field values are 
equal (B^ = -Be )■ A phase diagram for this situation 



With differentiating the equation of state (JSJ and @ 
with respect to the thermodynamic quantities T, B, n 
one can obtain generalized susceptibilities of the system, 
such as (dn/dfi)T,B or (dm/8B)T, n , etc. (the subscripts 
are referring to variables kept constant during differenti- 
ation). Focusing on the magnetic transition here we re- 
strict ourselves to the magnetic susceptibility of the sys- 
tem with fixed particle number, i.e. (dm/ dB)T, n - The 
region of interest in the (e s , B) plane is also restricted to 
the simply hatched region on Fig. and to the heavy 
line with (e,,. > e« , B = 0) , or equivalently to the sit- 
uations depicted on Figs. n a )~ c )i where one can find 
a purely magnetic transition above the BEC transition 
temperature. The temperature is also restricted above 
?BEc(e s ,-B). 

For this purpose let us first cast the expression of the 
density of the noncondensed particles © to a more ex- 
plicit form by performing the momentum integration: 



£(§) 

(2tt) 2 A 3 



F(^,f3[-n + c n n + r(mc s -gn B B)}), (11) 



where A = h/ ^j2Mk B T is the thermal wavelength of the 
atom, r(s) the gamma-function and F(s,j) is the stan- 
dard Bose-Einstein integral with parameter s and argu- 
ment 7 [23 ■ The derivative of i|ll|) with respect to B can 
be easily evaluated with the result: 

(12) 



~dB 



dm 
~dB 



where // = /1 — c n n was introduced for simpler nota- 
tion, and the relation dF(s, x)/dx = —F(s — l,x) was 
used |23|. Differentiating Eqs. © and with re- 
spect to B with n and T held fixed one arrives at a sys- 
tem of two equations for the quantities (dm/dB)x, n and 
(d/j,/dB)T.n-, from which the former one can be expressed 
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dm 
~dB 



with 



= 9Mb 



(P + R)Q + AP R 



P + Q + R + c s [(P + R) Q + 4P iJ] ' 

(13) 



P 



7? 



g£(|) 

(27r) 2 A 3 

flr(3) 

(2tt) 2 A 3 



F (§,/?[-// + (™c s - 9M bS)]), (14a) 
Ffo-P//), (14b) 

( o v , P(|,/3[-m' - (mc, - 5M bS)]) ■ (14c) 

The continuous magnetic transition is signalled by the 
divergence of the susceptibility Ijl3|) . or equivalently by 
the vanishing of its denominator: 

P + Q + R + c s [(P + R) Q + APR] = 0, (15) 

which is fulfilled along the thick curve of Fig. 
The magnetic susceptibility shows interesting properties 
also near the BEC. According to the expression Ijl3(l 
(dm/dB)T, n develops a cusp at T beo (-B) similarly to the 
susceptibility of the free Bose gas [13 at fixed number of 
particles. 



IV. CORRELATION FUNCTIONS OF 
GENERALIZED DENSITY OPERATORS 



Consider the following density operators 



^(k) = £(F 2 ) r , s „, 

q,r,s 



,t(k + q)a s (q). 



(16a) 
(16b) 



The former one is the particle density operator, while 
the latter one is the longitudinal magnetization density 
operator. Other generalized density operators can be de- 
fined as well for the spin-1 Bose gas, as done e.g. in Ref. 
pl|. but for the purposes of this paper the above two 
is enough. The correlation functions of the generalized 
density operators (|16f) are defined (for k ^ 0) as: 



(17a) 
(17b) 
(17c) 



with r being the imaginary time and T T the r ordering 
operator (see e.g. Ref. |2g). After the usual Fourier 
transformation into the Matsubara frequency represen- 
tation these correlation functions in general satisfy the 
following equations [2l| : 



Am(k,r) =-(T T [n(k,r)n t (k,0)] 
D zz (k, t) = -(t t [T z (k, r)^(k, 0)] 
D nz (k, t) = - h T [n(k, t)T\ (k, 



D nn — hYl nn 

d zz = nn zz - 



c n Tl nn D nri 

C n W nn D nz 



L nz &nz i 



c s Il nz D 



' zz i 



(18a) 
(18b) 
(18c) 



with the polarization functions: 

U nn (k,iuj n ) = y^II^(k, iu) n ), (19a) 

r,s 

n zz (k, lw „) = ^r S n^(k, 8Wll ), (i9b) 
n„ 2 (k,^„) = J2 S n^(k,iw„). (i9c) 

The polarization part 11^ (k, iui n ) is the contribution of 
interaction line irreducible Feynman diagrams, which can 
connect to an interaction line V^f from the right and to 
an interaction line V*'^ (a, b, c, d arbitrary spin indices) 
from the left |2l|. The interaction with four indices is 
defined in Eq. ©. (All correlation functions and po- 
larization functions in Eqs. H18|) depend on the (k, iu> n ) 
variables, which are omitted for the sake of brevity.) The 
solution of Eqs. I|18fl reads as: 



n _ fi n nn (i - c s ii zz ) + Cs u 2 nz 



det °s 



nz ~ det <V 



(20a) 
(20b) 



n _ fe Hzg(l c n Yi nn ) + c n Yi nz 

det % ' 1 UCj 

with the abbreviating notation 

det °| = (l - c„n„„) (i - Cs n zz ) - c nCs n 2 nz . (20d) 

The expressions H2()|l are quite general. They are valid 
both in the condensed and uncondensed phases and are 
contain no approximations. Of course the choice of the 
actual form of the polarization functions TF^ selects be- 
tween the different approximations and possible phases. 
The correlation functions are evaluated in the frame- 
work of the Random Phase Approximation (RPA) and 
for the uncondensed phase (n = 0), where the polariza- 
tion function is taken as the contribution of the bubble 
graph |2l| , which (for k and uj n not simultaneously zero) 
reads as: 



nn(k,tw n ) 



d 3 q 



'k+q,r 



(2n) 3 iuj n - h 1 (ek+q-ek) 



. (. 21 ) 

The expressions for finite magnetic field are quite similar 
to those of zero magnetic field |2l|. The effect of the 
magnetic field appears explicitly only in the contribution 
of the bubble graph l(2T|) and only through n q r , see Eq. 

©■ 



V. PROPERTIES OF WAVENUMBER 
DEPENDENT SUSCEPTIBILITIES 

The dynamics of the system can be studied with the 
help of elementary excitations. Spin dynamics is related 



6 



to longitudinal and transverse spin excitations. How- 
ever at the magnetic transition, where the susceptibility 
(dm/dB)T,n diverges longitudinal spin dynamics plays 
the relevant role as the soft mode of the transition. The 
spectrum of elementary excitations is related to linear re- 
sponse functions of the system (see e.g. ref. |2£|), which 
can be obtained from the correlation functions by ana- 
lytical continuation in frequency through the real axis 
from above. From now on we shall deal with such re- 
tarded correlation functions (response functions). With 
the help of linear response theory the density and mag- 
netization fluctuations of the system are related to the 
external potentials Sfi(k, ui) and 5B(k, ui) by 

Sn(k.,u) = — h~ 1 D nn (k, u>)S/i(k, ui) 

-g^Bh- 1 D nz {k,uj)SB{k 7 iu) 7 (22a) 
Sm(k, ui) — — h D zn (k, u))8fj.(k, ui) 

-g^Bfi- 1 D zz {k,u;)SB(k,Lu). (22b) 

From Eq. J25J: 

6(j,(k,w) = -g^B-^ — 7T7— T^ B ( k ' w ) ~ "■- 



D nn (k,u>) 



D nn (k,u) 



(23) 



Substituting it to Eq. (|22b|) one obtains the equation 



h 



5m(k, lu) 



H D nz {k,u>) 



Sn(k, ui) 



g^B 3MB D nn (k, ui) 

D nn (k, u))D zz (k,uj)~ D nz (k, ui)D zn (k,u>) 
D nn (k,u>) 



On this basis one can define 



6B(k,u). 

(24) 



_ ^_ 1 £>nn(k, u))D zz {k, lj) - D nz (k,uj)D zn (k,uj) 
D nn (k,u>) 

_ n nn (k,u})n zz (k,uj) ~ U nz (k,uj)Il zn (k,uj) 
II nn (k, uj)[1 - c s n zz (k, uj)] + c s II„ z (k, w)II z „(k, ui) ' 

(25) 

where Eqs. (|20[) were used to arrive at the last equation. 

The magnetic susceptibility ((13(1 can be expressed with 
the help of the response functions in the static limit (u> — 
0), since 



dm\ 6m(k,u>) 
= hm hm 



dBJ Tn Mw^o 5B(k,u) 
(with lim lim Sntk.ui) = 



(26a) 
(26b) 



In this limit the contribution to Sm comes from the sec- 
ond term of the r.h.s of Eq. 1)24(1 . According to the 
definition 1)25(1 it means that 



dm 



dB 



^) =-^limD£)(k,0). (27) 



T,n 



In the uncondensed phase the long wavelength limit of 
the static polarization functions take the following forms 
MM: 



Ilti ++ (k-0,0) 



n, 



n 



(r)- 



-P, 

-Q, 

(k^ 0,0) = -R, 



« 00 (k-0,0) 



(28a) 
(28b) 
(28c) 



and therefore IL m (k -> 0, 0) = -P - Q - R, Il nz (k -> 

0, 0) = -P + R and Ii zz {k -> 0,0) = -P - R. One can 
verify that the sum rule 1(27(1 is equivalent to Eq. ((13|) . 

When the density fluctuations are small the second 
term on the r.h.s. of Eq. ((241) will dominate for small 
k, ui values. (Note, that neglecting Sn, then S/j, and SB 
are related according to Eq. I(22a|) ~). Therefore we con- 
centrate on the expression 125|) of D zz . 

In the case when the momentum independent part of 
the Hartree energies (|6*cf can be assumed to be small, 

1. e. | — fi'\ -C ksT and \mc s ~ gnBB\ -C fcgT, and the 
wavenumber is much smaller than the thermal wavenum- 
ber feA < 1. Furthermore the frequency is also suffi- 
ciently small: fttku <C kX, the contribution of the bubble 
graph 1(21(1 can be approximated as pll l30l| 



w r r r (u,n) = 

+ in 



ft 



4tt 2 A 3 
1 



2,/oy 



Qu + i2J~a~ r 



with 



ay = — [3\p! + r(mc, 
u = kX, 

flu; fttku 
e k u z 



3(Ou + £2</oy) 3 
9HbB)], 



(29a) 



(29b) 
(29c) 

(29d) 



With these newly introduced notations the validity of the 
limiting form ((29a(l is equivalent to ay -C 1, |f2u| <C 1 and 
|n + i2yfor/u\ ^> 1. The expression l(29a(l is to be used 
in the formulas ((19(1 and ((2511 to obtain an explicit form 
for the response function D zz (k, u>). 

In the static limit one can set f2 = in the expression 
5a| safely with the result: 



n£(u,o) 



ft 



4"7T 2 A 3 



24a 



3/2 



(30) 



which can be used in the polarization functions ((19(1 to 
obtain 



n, 



n 2 



-P-Q-R; P = P-pu 2 



ft 



1 



96ttA 3 ct 3 / 2 ' 



-P + R: 



-P-R; 



Q = Q — qu 



R = R-ru 2 



ft 



1 



96ttA 3 3 / 2 ' 



1 



96ttA 3 ct 3 / 2 
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12 1.421 1.422 1.423 1.424 1.425 1.426 1.427 

T c /To T/T 

FIG. 3: The damping rate —V of the longitudinal spin ex- 
citations (spin-density fluctuations) as a function of T/To at 
B — B c and e s = 0.9 (upper panel). The damping rate is 
given in units of ksTo/h. The magnetization density is plot- 
ted on the bottom panel for enlightenment. 



up to the order of u 2 . Therefore the long wavelength 
limit of the static susceptibility at T — T c , where Eq. 
(|15|l holds takes the following form: 



i>W(fc )£ j = o) = 



(fcA) 2 

Q{P + R) + APR 



p + q + r + c s [q(P + R) + Q(p + r) + A(Pr + pR)] ' 

(31) 

From Eq. l)."ilH one can see that the critical exponent 77, 
defined as C(k) oc k~ 2+n at the critical temperature for 
the correlation function C of the order parameter, is zero 
(77 = 0). 

To obtain the spectrum of the longitudinal spin fluc- 
tuations one has to look for the poles of the response 
function l|25[) . or equivalently look for the zeroes of its 
denominator. The resulting equation is as follows: 

Il„ n (k, to) [1 - c s II zz (k, u)] + c s Tl nz (k, oj)U zn (k, uj) = 0. 

(32) 

With the help of the numerical approximation of the an- 
alytic continuation of the bubble graph (|21[1 30l l3l| and 
with the help of Eqs. ((TH . Eq. (|3^j) can be solved di- 
rectly. In the long wavelength limit one finds a purely 
imaginary (overdamped) excitation spectrum depending 
linearly on wavenumber: 



u(k) = —iT k, 



(33) 



with r the velocity dimensional damping rate. The nu- 
meric values of T are given in Fig. [3]for e s — 0.9, B = B c 
and temperature values above that of the BEC. One can 
clearly see the soft mode (r = 0) at B = B c and T = T c . 

For analytical calculation at B c and near T c one can 
suppose that \flu\ <C 2^/5^ for all r, since the mode is 
soft and neither of the as is zero. In this case Eq. <|29a[) 



can be further approximated. Taking only the first three 
terms form the expression and casting the last one into 
series of flu/2^/dr one arrives to: 



E£(u,n) = - 



ft 



4tt 2 A 3 



-F(±,<7 r ) +i 



flu 
4ov 



(34) 



Using Eq. I)29d(l and for T one can obtain 

r = A 

0H 

P + Q + R + c s [Q{P + R) + APR] 
p + q + r + c s [q{P + R) + Q(p + f) + A{P? + pR)] ' 



where P, Q, R are given by Eqs. 114|) and 

ft 1 



P 



167rA 3 fi+ ' 
ft 1 

167rA 3 <7q ' 

_l 1_ 

167rA 3 (7- 



(35) 

(36a) 
(36b) 
(36c) 



The numerator of the damping parameter (|35|l is equal 
to the denominator of the static, homogeneous suscepti- 
bility l)13|l . Accordingly the damping rate, T, vanishes at 
T = T c (for B = B Cl where the transition is continuous) 
and r is proportional to the inverse of the susceptibility 
(|13|l for T T c , like in conventional theory. 



VI. CONCLUSIONS AND SUMMARY 

We have studied the statistical physics of the spin-1 
Bose gas with ferromagnetic coupling in the presence of a 
magnetic field both from the static and dynamical points 
of view. The equation of state and the phase diagram of 
the system was investigated both for the Bose-Einstein 
condensed and uncondensed phases. The discussion of 
other static quantities and the dynamics is concentrated 
to the noncondensed phase, where the purely magnetic 
phase transition takes place, which occurs prior the Bose- 
Einstein condensation if the spin-flipping coupling con- 
stant (|c s |) is large enough. 

This magnetic transition was found to be of first or- 
der in general. There is always a critical magnetic field 
value, where the magnetic transition is continuous and 
above that there is no phase transition at all (magneti- 
zation behaves analytically). There is also a tricritical 
point (at e s = 6s ~ 1.22, T w 1.58), when the criti- 
cal magnetic field reaches zero. Above this value of e s 
the usual paramagnetic-ferromagnetic phase transition 
occurs in zero external magnetic field. The static sus- 
ceptibility (dm/dB)T,n shows a divergence at the critical 
points of the continuous MT, while it has only a cusp at 
the point of the BEC. 
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The long wavelength dynamics also exhibits features 
of the MT. The longitudinal spin excitation becomes a 
soft mode in accordance with the diverging susceptibility 
in the vicinity of the critical point of the MT. Moreover 
the soft mode behaves like in conventional theory; close 
to T c the damping rate is proportional to the inverse of 
the susceptibility. 
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